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PREFACE 


This report is expected to receive wider distribution than any 
published previously on the same study. It therefore incorporates much 
of the material contained in those prior reports, for the convenience of 
readers who have not had access to them. Those already familiar with 
the study concept and prior progress are entreated to understand and 
forgive the duplications. 

The study is the first part (a feasibility investigation) of a longer- 
range application development project whose concept has evolved from an 
idea suggested at the International Astronautical Congress in Brussels in 
1971 . Since it was initiated under NASA contract in July 1971, it has 
focused on the same objective even though (as in all such research) the 
technical scope, methodology and manpower levels have from time to time 
been revised. A brief history of the study appears in Appendix A. 

The study objective and concept are described in Section 1, with a 
summary of results. Section 2 supports the choice of simulation model 
and describes its use in the study. Section 3 presents the study methodology 
and describes each task and the results of each task. Recommendations for 
additional research and validation efforts appear in Section 4. Appendices 
contain a summary project history, a list of references and the basic 
distribution list for the report. 

The authors are grateful to many individuals for initiation and support 
of and contributions to the study. Dr. Peter A. Castruccio (presently 
president of Ecosystems International, Inc. ) and Mr. Andrew Adelman of IBM 
were energetic in promoting the basic idea. Dr. George McDonough of 
NASA's George C. Marshall Space Flight Center sponsored the project 
and was instrumental in securing funds to support it. Dr. William D. Clarke, 
as NASA-MSFC technical manager, made several helpful suggestions and 
first recognized the potential value to others of the processed and integrated 
data base which is a by-product of the study. 

The advice of authorities in many areas related to streamflow fore- 
casting and hydrologic modeling was available to the study team. 

Dr. L. Douglas James of the Environmental Research Center, Georgia 
Institute of Technology, provided the basic simulation and calibration 
programs, advised us in their implementation and reviewed interim results 
of the study. In obtaining historical and physiographic data, selecting 



watersheds for study, and understanding the objectives and methods of 
operational river forecasting, the authors have received invaluable 
assistance from personnel of the Tennessee Valley Authority; National 
Weather Service Hydrology Office and Lower Mississippi River Forecast 
Center; U. S. Geological Survey, in several locations; U. S. Soil Conser- 
vation Service; and the U. S. Army Corps of Engineers. 

Mrs. Becky Robinson of IBM deserves a special measure of thanks. 
In a remarkable exhibition of versatility, in addition to performing as 
departmental secretary, she: 

• secured storage facilities for the tons of data acquired 
and generated in the study; 

• constructed photomosaics and overlay graphics; 

• identified and maintained hundreds of maps and aerial 
photographs; 

• made special trips to TVA in Knoxville to obtain 
historical data; 

• prepared illustrations for presentations and reports; 

• typed, edited and assembled technical reports, including 
this one. 


These acknowledgements of assistance do not imply any endorsement 
or approval by those whose help we appreciate so much. The conduct of 
the study, the results achieved and the content of this report are the 
responsibility of the authors. Any errors of fact, methodology or logic 
should be attributed solely to us. 




Reuben Ambaruch 
Principal Investigator 





John W. Simmons 
Study Manager 
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SECTION 1 


INTRODUCTION 


1.1 OBJECTIVE 


The objective of the study whose results are reported herein is to 
assess the feasibility of using the data produced by remote observation from 
space and/or aircraft to reduce the time and expense normally involved in 

achieving the ability to predict the hydrological behavior of an ungaged 
watershed. Such a capability will enhance effective planning for urban and 
industrial development, flood control, hydroelectric power, navigation, 
and water resources management. Traditionally, the ability to predict 
streamflow rates, in response to given precipitation, is attainable only after 
installation of extensive instrumentation and collection of data for several 
years to estimate the operating parameters of a model of the basin or water- 
shed system. This project was aimed at developing techniques for using the 
large amounts of data acquired by earth observation missions to determine 
those parameters more directly. 

1.2 BASIN MODELING AND STREAMFLOW FORECASTING 

That aspect of hydrology known as Streamflow Forecasting undertakes 
to predict the outflow from a given catchment, in terms of flow rate as a 
function of time, in response to a given precipitation event under given 
initial conditions. This capability is vital to effective planning for urban/ 
industrial development, flood control hydroelectric power, navigation, and 
water resources management. 

Figure 1 depicts the cross section of a somewhat idealized rural 
catchment and identifies the principal phenomena at work in the rainfall- 
runoff relationship. The input (precipitation) is partially intercepted by 
vegetation and water retention areas. Moisture reaching pervious surfaces 
divides between overland flow, infiltration and evaporation. Through sub- 
surface processes, interflow and groundwater flow contribute ultimately to 
streamflow, with some losses due to transpiration through plant life. 

All the phenomena involved in the hydrological cycle are widely and 
well understood qualitatively, and several empirical relationships have 
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Figure 1. Cross Section of Idealized Rural Catchment 



been developed. Integrating them into comprehensive parametric models 
has been achieved by several investigators, developments which were 
facilitated by the availability of large, high-speed digital computers. Of 
these, the Stanford Watershed Model is probably the best known. In order 
to apply any model to a given catchment, however, it is necessary that it 
be calibrated . That is, several years of streamflow are simulated using 
actual precipitation data, synthesized flows are then compared to actual 
recorded flows, model parameters (as well as seasonal factors) adjusted 
by hydrologists and the process repeated until acceptable simulation 
accuracy has been achieved. More recently, self- calibrating models have 
been developed which use trial and error routines to optimize model and 
seasonal parameters, but the calibration still required historical stream- 
flow data. If a stochastic model or technique is used instead of a parametric 
model, several dozen storm events in the form of previously measured rain- 
fall and runoff data must be used in multivariate statistical analyses in order 
to establish confidence in the results. 

Whatever the model or analytical techniques used, current methods 
of streamflow forecasting require installation and maintenance of a system 
of instruments so that data may be collected for several years. The basin 
of a major river must be divided into hundreds of watersheds of manageable 
size and each one provided its own stream gage. Although the installation 
and maintenance expenses may not be objectionable, waiting for several 
years of data collection to initiate a new development may be intolerable. 
There is a real need for a method of predicting the hydrological behavior 
of a previously ungaged catchment quickly (within two months, say, rather 
than three to five years) and with reasonable accuracy. 

1.3 POTENTIAL APPLICATION OF REMOTE SENSING 


There has been increasing interest in recent years in the relationship 
between hydrological behavior of a watershed and a unique set of physio- 
graphic features by which it might be characterized. One would like to be 
able to quantify simulation model parameters, either directly or inferentially, 
from observations of catchment morphology and climate, taking advantage of 
all available prior knowledge of the area. Earth observation spacecraft, 
manned or unmanned, make observations using a variety of sensors, cover- 
ing large areas in a short time. The objective of the study whose results 
are reported herein is to assess the extent to which the data provided by 
earth observation missions can be used to reduce the time, effort, and 
cost of predicting the hydrological behavior of a given catchment. 
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1 . 4 FEASIBILITY STUDY METHODOLOGY 


The study was conducted in three phases as shown in 
Figure 2, aimed at devising a set of widely applicable pre- 
diction model parameters in terms of photographically observable 
characteristics (such as climatology, areas, elevations, and land use) 
and inferable characteristics (such as soil depth and porosity). 

1.4.1 PRINCIPAL STUDY ELEMENTS 

The principal elements of the study are a simulation and a cali- 
bration model, a study area, and a complete historical and physiographic 
data base pertaining to the study area. They are as follows: 

• Simulation Model. The continuous simulation model chosen 
for the study, for the reasons discussed in Section 2, is 
known as the Kentucky Watershed Model (KWM). This model 
is based upon the well-known Stanford Watershed Model IV, 
and has been further adapted and refined by IBM for application 
to this project. 

• Calibration Program. A calibration program known as OPSET 
(because it estimates the optimum set of watershed model 
parameters) is a companion to the KWM. It has also been 
refined and modified for application to this study. By atrial 
and error process it estimates values for 13 simulation 
model parameters which are further refined by manual 
operations. This program is also described in Section 2. 

• Study Area. The area from which test watersheds are chosen 
for the study is the Tennessee River Valley, a major watershed 
of approximately 104,000 square kilometers in the southeastern 
United States. This is a well-instrumented, thoroughly photo- 
graphed and mapped area for which copious historical data records 
are available. 
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Figure 2. Study Methodology 










Within the Tennessee Valley 55 watersheds have been identified 
as suitable to the purposes of the project, representing all six 
physiographic provinces of the valley. It was estimated that 35 
watersheds would be needed for confidence in the results of the 
correlation process. Of these, 25 were to be used for calibration 
and 10 for tests of simulation accuracy. Table 1 lists the 
limited set actually used in the study. Figure 3 indicates 
their locations in the region. These represent five of 
the six physiographic provinces and range in area from 7 to 365 
square kilometers. The average overland flow surface slopes 
are indicative of the varying ruggedness of the region, from 
the Nashville Basin to the Blue Ridge Mountains. For all except 
the smallest watersheds, hourly rainfall station records are 
supplemented by daily rainfall records for the sake of accuracy. 

• Historical Data Base. The study requires accumulation of a 
comprehensive data base consisting of streamflow, precipitation, 
and evaporation records for the study area. The Tennessee 
Valley watershed is one which is extensively instrumented and 

for which ample historical data are available. It contains approxi- 
mately 560 rainfall stations, both hourly and daily. Almost 1,000 
streamflow records have been accumulated in the past several 
decades; approximately 360 streamflow stations are currently 
active, and some 600 have been discontinued. Other climato- 
logical data such as temperature, evaporation, wind and humidity 
have been carefully collected and catalogued. 

• Physiographic Data Base . The study also requires the accumulation 
of a comprehensive physiographic data base for the study area. 
Complete topographic map coverage is available, consisting 
almost entirely of 7. 5 minute quadrangle maps derived from 
aerial photographs at a scale of 1:24,000. The Tennessee 

Valley has also been completely surveyed by aerial photography, 
and black and white photographs are available for every water- 
shed, mainly at a scale of 1:24,000. Some photographic coverage 
is also recently available in pseudocolor infrared taken from 
high altitude aircraft, at a scale of 1:63,000. Excellent images 
from the Earth Resources Technology Satellite (ERTS) have also 
been obtained. 


1-6 



1-7 


Table 1 . Selected Watersheds 



CATCHMENT NAME 

PHYSIOGRAPHIC 

AREA 

(SQUARE 

AVERAGE 


u # 

s 



PROVINCE 

(KILOMETERS) 

SLOPE 


Daily 

E 

1 . 

WHITE CREEK NEAR SHARPS CHAPEL, TN (WHITE HOLLOW) 

Valley and Ridge 

6.94 

.300 

1 

- 

C 

z 

LITTLE CHESTUEE CREEK BELOW WILSON STATION, TN 

Valley and Ridge 

21.30 

.237 

1 

1 

C 

3. 

SOUTH FORK MILLS RIVER AT PINK BEDS, NC 

Blue Ridge 

25.70 

.315 

2 

1 

C 

4. 

NOLAND CREEK NEAR BRYSON CITY, NC 

Blue Ridge 

35.70 

,500 

2 

- 

V 

5. 

WEST FORK PIGEON RIVER ABOVE HAZELWOOD, NC 

Blue Ridge 

37.30 

.413 

1 

1 

C 

6. 

BIG BIGBY CREEK AT SANDY HOOK, TN 

Highland Rim 

45.30 

.185 

1 

1 

V 

7. 

BIG ROCK CREEK AT LEWISBURG, TN 

Nashville Basin 

64.50 

,155 


- 

C 

a 

LITTLE RIVER ABOVE HIGH FALLS, NEAR CEDAR MTN., NC 

Blue Ridge 

69.50 

.246 

1 

2 

V 

9. 

TRACE CREEK NEAR DENVER, TN 

Highland Rim 

78.70 

.129 

1 

2 

c 

10. 

PINEY CREEK NEAR ATHENS, AL 

Highland Rim 

145.00 

.033 

1 

4 

V 

11. 

RUTHERFORD CREEK NEAR CARTERS CREEK, TN 

Nashville Basin 

178.00 

.133 


2 

u 

1Z 

EMORY RIVER NEAR WARTBURG, TN 

Cumberland Plateau 

215.00 

.298 


3 

c 

13. 

SEWEE CREEK NEAR DECATUR, TN 

Valley and Ridge 

303.00 

.165 ! 

2 

2 

C 

14. 

TOWN CREEK NEAR GERALDINE, AL 

Cumberland Plateau 

365.00 

.062 

2 

5 

C 

15. 

POPLAR CREEK NEAR OAK RIDGE. TN 

Cumberland Plateau 

214.00 

.329 

1 

3 

c 

16. 

EAST FORK POPLAR NEAR OAK RIDGE 

Cumberland Plateau 

50.50 

.215 

1 

1 

V 


USE: C • CALIBRATION; V = VALIDATION; U = UNUSE ABLE (INADEQUATE DATA) 
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1.4.2 CALIBRATION 


In the calibration phase the model parameters and their related 
observables were quantified and adjusted for ten selected watersheds 
from five of the six physiographic provinces of the Tennessee Valley. 
These sets of satisfactory model parameter values are those which 
yield simulated streamflows which best correlate with actual recorded 
flows for periods of five to ten years. A sensitivity analysis of two 
watersheds provided guidance for final manual adjustments of the model 
parameters, thereby minimizing the time and effort spent in the cali- 
bration phase. 

1.4.3 CORRELATION 

In the correlation phase, satisfactory model parameters were 
correlated with observable physiographic characteristics in a multiple 
regression analysis, using the results of ten calibrations as a data base. 
The product of this phase is a set of regression equations, each of which 
is an expression for a model parameter as a linear combination of other 
parameters. 

1.4.4 VALIDATION 

In the validation phase, the relationships produced by the correlation 
analysis were tested on five watersheds from different physiographic 
provinces of the region. No calibrations were performed on these basins; 
they were treated as ’’ungaged” watersheds. Their simulation model 
parameters were quantified by using the regression equations, and then 
the simulation model was run, using those parameters and historical 
climatological data as input. The simulated streamflows compared well 
with observed streamflows, in respect to most indices of performance 
indicating the feasibility of the technique. 

1. 5 SUMMARY OF RESULTS 

1. 5. 1 FEASIBILITY CLEARLY INDICATED 

The five validation runs produced simulated streamflows which 
correlated remarkably well with observed streamflow. Daily correlation 
coefficients ranged from 0. 83 to 0. 87; monthly, from 0. 92 to 0. 97. Many 
major storms were reasonably well matched wi th respect to peak flows 
and timing of peaks. For a multi-year open-loop simulation, this is 
adequate for most applications, and it strongly indicates the feasibility 
of using remotely sensed data to forecast the hydrologic performance of 
an ungaged watershed. 
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This happy conclusion should be kept in the perspective shaped by 
the following facts and circumstances. 

• The initial study plan envisioned using 25 watersheds for 
calibration and ten for validation to provide statistical 
confidence in the results. The study actually used ten 
and five respectively, leaving some (undefined) un- 
certainty in the validity of the result. 

• The remote sensing imagery used to determine watershed 
physical characteristics was actually aerial photographs, 
mostly at a scale of 1:24, 000. No interpretation of 
higher -altitude (RB-57F, U-2) or space (ERTS) images 
was done in the study. 

• Some ground truth was necessary, in this instance soil data. 

This is not considered a disqualifying fact, because the need 
for some ground truth probably cannot be obviated in future 
applications. 

• The correlation (i. e. , multiple regression) analysis deserves 
some refinement, based on an enlarged appreciation by the 
statistical analyst of the hydrologic processes and relationships 
found in a watershed. 

1. 5. 2 SIGNIFICANT BY-PRODUCTS 

There are two by-products of the study which the authors believe 
to be of potential value to other investigators in hydrologic modeling and 
related fields. 

1. 5. 2. 1 Integrated Data Base 

As in most applications using simulation models, the preponderant 
share of the study effort was devoted to acquiring, verifying, formatting, 
and preprocessing input data to construct a data basis. As a result of this 
study, a data base now exists in the form of a magnetic tape. The tape has 
a capacity of up to 50 files, one file per basin. Each file has fields allocated 
for up to ten years of historical climatological data (mean basin hourly 
rainfall, evaporation, temperature), mean daily streamflow, and selected 
storm events (up to ten per year, characterized by peak flow and time of 
peak). Each file also has fields for watershed features, simulation model 
parameters and program control options. 

At present, the data base contains data, except temperature, pertinent 
to the 15 Tennessee Valley basins used in the study. To any other investigator 
interested in doing similar studies in the same area, the integrated data base 
may be of great value. 


1-10 



1. 5. 2. 2 The IBM System for Simulation and Analysis of Watersheds 


Since the study team began adapting the Kentucky Watershed Model 
and its companion calibration program, OPSET, to the feasibility study, 
several programming additions and modifications have been made. The 
result is a system of methods and software, built around the basic simu- 
lation and calibration models, which is a powerful hydrologic research tool 
readily adaptable to operational use. Features of the system are as follows. 

• Input data preprocessing and verification. 

• Plot outputs. 

• Tabular summaries. 

• Storm analysis. 

• Statistical analysis routines. 

• Terminal operation for 

set up and initiation of calibration/simulation runs, 
with a batch-processing computer, or 

real-time interaction, with a time sharing computer. 

The effort required to develop this system for efficient operation has 
been more than compensated by the manpower saved in manipulation of 
data, handling many large decks of cards, and analyzing raw simulation 
outputs. In fact, the study could not have been completed without it. 

1.6 A GLIMPSE OF THE FUTURE 

Section 4 of this report summarizes research efforts worthy of con- 
sideration as sequels to this feasibility study. Several years will be required 
for the idea tested preliminarily here to be reduced to practice in another 
region. Nevertheless, one can readily visualize the sequence of events. 

• A set of spacecraft sensors observe an ungaged region con- 
current with aircraft underflights and ground truth patrols. 

• The data so acquired are analyzed and interpreted to yield 
watershed characteristics. 

• Statistically derived relations are used to quantify simulation 
model parameters from observed characteristics. 
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The simulation model is run, using climatological inputs 
typical of the region, to generate multiyear streamflow 
predictions and statistics, for planning purposes. 




This can be done in a short time, two to three months, without extensive 
instrumentation or the need to take calibration data for three or more years. 
Residual benefits will derive from implementing the system for periodic 
short-term river forecasting, using a few well-placed stream gages for 
"closed- loop” operation and a forecasting accuracy within five per cent. 
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SECTION 2 


SELECTION OF SIMULATION MODEL 


2.1 OBJECTIVE 


The spreading availability of high-performance computers has 
stimulated remarkable progress in development of hydrological models. 
An individual investigator may prefer one model or another on the basis 
of several criteria, with emphasis on the objectives of his research. We 
have selected a model with which we can accomplish the calibration, 
correlation, and validation phases within the study schedule time and 
manpower budgets. The models and procedures developed must be 
capable of predicting the hydrological performance of an ungaged basin 
with parameters derived from space or aerial imagery, sparse ground 
samples, precipitation and evaporation data. 

2.2 STOCHASTIC MODELS 


The project first considered stochastic models, which involved the 
use of multivariate regression analysis to develop predictions of runoff 
as a function of a limited number of observable variables: storm duration, 
storm intensity, time of the year, and initial moisture conditions — the last 
estimated empirically from recent precipitation experience. Techniques 
were identified for deriving rainfall- runoff relations using observable or 
inferable physiographic features. Synthetic hydrograph* methods were 
considered for determining streamflow as a function of time. Methods 
were identified for statistically correlating observable physiographic 
features with unit hydrograph parameters. 

2.3 PARAMETRIC MODELS 


Parametric models were also considered and found to be preferable 
to stochastic models for the following reasons: 


* 

A hydrograph is simply a plot of streamflow in volume per unit time or 
river height as a function of time. See Reference 2, Chapter 9. 
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1 . The cause and effect relationships among the conditions and 
processes in a watershed are obscured in the stochastic model 
but are used explicitly in the design of the parametric model, 
many of whose parameters are observable directly from remote 
sensed imagery. 

2. In using a stochastic rainfall- runoff (RF-RO) model, it is 
necessary to analyze some 50 observations (storm events) per 
watershed to obtain confidence in the validity of the derived 
relationship. In determining RF-RO relationships based upon 
physiographic parameters, approximately 50 watersheds are 
needed, each one having 50 associated storm events. 

3. In order to take into account all practical variations in and 
interrelations among watershed parameters, a prohibitive 
amount of data is required for thorough statistical analysis. 
Thomas and Benson** identified 71 physiographic and statistical 
variables of interest in relating streamflow characteristics to 
basin physical characteristics. 

4. A parametric model that is capable of simulating continuous 
flows over a long period of time supports a larger variety of 
planning objectives than a method which is intended to estimate 
or predict a single variable or isolated event (e.g. , flood peak 
or total runoff from a storm). 

2.4 STANFORD WATERSHED MODEL (SWM) 

4 5 

The Stanford Watershed Model ’ is probably the best known of the 
parametric hydrological models and, in all its modifications, is probably 
the most widely used. Since it was originally published in 1962, several 
reports have appeared in literature describing modified versions and 
applications (References 5, 6, 7, 8, 9, 10, and others). As a proven 
tool it was attractive to our project, whose scope does not embrace 
advances in hydrological modeling. 


The Stanford Watershed Model uses a moisture accounting system 
to synthesize a continuous hydrograph from the following: 
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1. Recorded climatological data, precipitation, evaporation and 
(for snowmelt situations) temperature, 

2. Measurable watershed characteristics such as drainage area 
and friction of the watershed in impervious surfaces, and 

3. Parameters used in the computation process which are known 
to vary in magnitude among watersheds but have not been 
quantitatively tied to specific measurable watershed properties. 
For example, one parameter indexes the capacity of the soil 

of the watershed as a whole to retain water. 

The third class of inputs requires a trial and error series of cali- 
bration runs to quantify a set of model parameters which will synthesize 
flows with acceptable accuracy. 

Figure 4 depicts the accounting of moisture entering the watershed 
until it leaves by streamflow, evapotranspiration, or subsurface outflow. 

A series of relations each based on empirical observation or theoretical 
description of a specific hydrologic process, is used to estimate rates 
and volumes of moisture movement from one storage category to another, 
in accordance with current storage states and the calibrated watershed 
parameters. The model routes channel inflow from the point where it 
enters a tributary channel to the downstream point for which a hydro- 
graph is required. A subroutine exists for including snow in the accounting 
but is not needed in the present study. 

2.5 KENTUCKY WATERSHED MODEL (KWM) AND OPSET PROGRAM 

The Stanford Watershed Model was originally written in the Burroughs 
Computer Language (BALGOL) then in use at the Stanford Computer Center. 
It has subsequently been translated into Fortran IV, and a number of 
adaptations were introduced in one version to suit the climate and 
geography of Kentucky as representative of the eastern United States. 

In a recent research program a version of the model using an initial 
set of model parameter values and a number of control options was 
developed for use with a self- calibrating streamlined version of the 
model. These models are referred to as the Kentucky Watershed Model 
(KWM) and OPSET (because it estimates the Optimum SET of model 
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Figure 4. Moisture Accounting in the Stanford Watershed Model 
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parameters). The availability and utility of these models and reports 
describing them 13 [ e£ j to their use in the IBM project. 

2.5.1 KENTUCKY WATERSHED MODEL 

Figure 5 lists the inputs (exclusive of control options) used by the 
Kentucky Watershed Model to simulate streamflow. Climatological data 
can be obtained from precipitation records or can be hypothetical, the 
latter being useful in generating rainfall-runoff predictions. The inputs 
classed as "Overland Flow Parameters" and "Watershed Parameters" 
are readily obtainable from analysis and interpretation of images (maps 
and/or photographs). The inputs on the right side of the figure are 
estimated in the calibration phase by OPSET. Some additional manual 
calibration is necessary to develop a set of model parameters that best 
represents the watershed. 

2.5.2 OPSET 

When a user applies a simulation model to a watershed, there are 
several parameters whose values he must initially guess and subsequently 
adjust, between trial runs of the model and comparisons of synthesized 
with observed flows. This trial and error calibration r equires ingenuity, 
familiarity with how the parameters interact within the model, and some 
understanding of the sensitivity of simulated flows to specific parameter 
adjustments. The process is aided greatly by a thorough understanding 
of the hydrologic process and by the guidance published by Crawford and 
Linsley^**. Through careful parameter adjustment, one can cause sim- 
ulated flows to approximate recording flows but never to match them 
exactly. Several combinations of parameter values can produce com- 
parable results from an overall viewpoint, and the final choice may well 
hinge on whether a particular comparison emphasizes flood peaks, annual 
runoff volume, or some other hydrograph feature. The final acceptance 
of a set of parameters may depend heavily on subjective factors. 

12 

In developing OPSET, Liou provided a tool for calibrating the 
KWM with a minimum of subjective decisions. The parameter optimi- 
zation concept is depicted in flow chart form in Figure 6. The input data 
consists of control options and initial conditions as well as the inputs 
listed in Figure 7. A simulation is performed, one year at a time, using 
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Figure 7. Calibration Program (OPSET) Inputs and Outputs 




a "streamlined" KWM. The synthesized flows are compared with the 
observed flows. An objective function is used to determine when an 
optimum set of parameters has been found. If the best match has not 
been achieved, parameters are again adjusted and the simulation run 
again. This sequence is repeated until a satisfactory parameter set has 
been quantified. 

Figure 7 also lists the 13 outputs of OPSET, in addition to simulated 
streamflow. (Comparison with Figure 5 shows the relationship to KWM. ) 
These parameters are the most difficult to measure directly and ones to 
which simulated flow values are sensitive. Yet they, too, must be esti- 
mated from observation and correlation in the future, to apply the 
simulation model to an ungaged watershed. The calibration process 
should be based on three separate water years for the same basin. 
Simulation model parameters are then derived by averaging the results 
of the three calibration runs. A minor modification to OPSET has been 
implemented to generate a more precise Base Flow Recession Constant 
(BFRC). As it is presently designed, OPSET estimates parameters 
which produce accurate simulations of major winter storms (with 
respect to flood peak magnitude and timing) but misses summer and 
autumn storm peaks by significant factors. Manual adjustments are 
required to achieve accurate simulation in the latter. An improvement 
in OPSET efficiency could be achieved by modifying it to calibrate on 
the basis of several consecutive years rather than one year at a time. 

2.6 MODIFICATIONS AND ADDITIONS TO KWM AND OPSET 


In adapting the KWM and OPSET programs to the project, several 
additions and modifications were implemented by IBM for improved 
flexibility and efficiency, both in running the program and evaluating 
results. They are listed as follows: 

1 . The subroutine called READ was modified to operate on any 
current System/360 installation. 

2. Both models are in the Production Library, and the input data 
is in a partition data set (PDS) based operation which will (a) 
avoid the need for maintaining several thousand cards for each 
watershed, and (b) permit operation from a terminal, giving the 
investigator quick access to the data base for modifications. 
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3 . 


A subroutine designated INT has been added to convert 
historical stream gage height information from strip charts 
to streamflow rates (volume per unit time), applying the 
gage rating tables. It also analyzes the peaks and calculates 
the runoff. 

4. A subroutine designated STAT has been added to provide 
statistical analysis of monthly and daily data for a single 
year, and multiple years. Sample output for a yearly 
statistical summary is shown in Figure 8, and the total 
statistical summary output for seven years of simulation 
is presented in Figure 9. 

5. An Output Summary routine tabulates all the end results of 
a calibration or simulation run for a particular watershed 
so that the investigator need not search through several 
pages of printout to compare and evaluate results. Figure 10 
presents an Output Summary. 

6. As it was originally received, the KWM could perform 
simulations only for one year at a time, and it was necessary 
to re-establish initial conditions as well as re-read watershed 
descriptive data before each run. The capability has been 
introduced to simulate any number of successive years with 
only one reading of the watershed description and retaining 
simulated moisture conditions at the end of each year to 
serve as initial conditions for the succeeding year, thereby 
providing more accurate simulation results. 

7. A Plot Output routine with several options was introduced in 
order to minimize manpower and time needed for evaluation 
of results, sensitivity analysis and calibration "fine tuning." 
Plots shown in a subsequent section of this report were pro- 
duced by this program. Vertical (flow rate) and horizontal 
(time) scales may be specified or automatically chosen. Plots 
are available for simulated flow, observed flow or the two 
superimposed on the same scale. Time intervals available 
are one year (mean daily streamflows), one month (mean daily 
streamflows) and storm duration (hourly streamflows for the 
duration of the selected storm runoff). 
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The subroutine INT, added to KWM for storm analysis, constitutes 
a major modification. It has the following capabilities: 

• Converts gage height (ft) to streamflow (cfs), a task which 
formerly had to be done manually. 

• Calculates the total runoff volume resulting from each storm. 

• Summarizes the characteristics of each selected storm: peak 
flow rate, time of peak, etc. 

Evaluation and analysis of hydrograph characteristics requires 
observed hourly streamflow data. An example of available observed* 
hydrograph data is shown in Figure 11. The observed hydrograph data 
is in terms of stream gage height (ft), but the simulated hydrograph data 
is in streamflow (cfs). Therefore the observed data must be converted 
to cfs, which requires rating tables for for that period of time. A sample 
rating table is shown in Figure 12. This rating table is valid for the period 
from February 11, 1953, to March 9, 1953. The present study uses seven years 
of historical data (1950-1956) and requires 17 rating tables. The subroutine 
INT uses the proper rating table to convert stream height to flow rate. 

Operating with selected major storms, rather than a continuous 8,760 hours 
of data per year, reduced the effort associated with observed input data 
(manual translation of gage height hydrographs) and processing time by a 
factor of 40, assuming an average storm duration of 48 hours and five 
major storms per year. The subroutine INT has the capability of integration 
to calculate the runoff volume produced by observed and simulated storms. 

The calibration and simulation models and their data base have been 
adapted to operation from data terminals in the IBM- Hunts ville facility. 

Access to an IBM System/360 Model 75 is afforded from a Model 2260 
terminal, from which an operator can call up sections of the data base 


^Throughout this report, the term ’’observed" refers to recorded historical 
data, distinguished from the "simulated" or "synthesized" results produced 
by the simulation model. 


2-14 



*or»7A L |.: ,*AiLfl avu ;:CvBO«i» 



Figure 11. Example of Observed Hydrograph Data 
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Figure 12 . Typical Stream Gage Rating Table 
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and modify elements thereof, enter control options in the program, 
change model parameters, and enter the program in the computer 
system queue, without having to manipulate punched cards. In the 
sensitivity analysis (explained in a subsequent section of this report), 
several runs can be set up and initiated at one setting. It is conserva- 
tively estimated that this terminal operation approach has obviated 
the storage and maintenance of 500,000 punched cards and improved 
the efficiency of operation by a factor of ten. 
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SECTION 3 


TECHNICAL RESULTS 


This section describes the principal and secondary accomplishments 
of the feasibility study. The study consists of nine tasks, identified as 
shown in Figure 13, which illustrates the flow of study activities. In the 
remaining pages of this section, each task is described with respect to 
its objective, the technical approach used, and the results achieved. 

Physiographic data were collected (Task 1) in the form of topographic 
maps, aerial photographs, and other documents for 15 watersheds of the 
Tennessee Valley. The historical data (Task 2) consists of observed 
streamflow, precipitation, and evaporation data related to the same 1 5 
watersheds. All these physiographic and historical data are digitized 
and integrated into the Master Data Bank (Task 3). (Note: Although only 
15 basins were used in the study, the study team’s accumulation of raw 
data pertains to some 50 basins. ) 

Ten of the test watersheds were modeled, and the models were 
calibrated using OPSET and manual adjustments (Task 4) in order to 
derive a set of optimum model parameters for each watershed. This 
process was aided by the sensitivity analysis (Task 5), which was per- 
formed on two of the selected watersheds. The sensitivity analysis 
shows quantitatively how the accuracy of the simulation model is affected 
by variations in the model parameters. 

The physiographic data base yields watershed characteristics for 
each of the 10 basins on which a calibration will have been performed. 

These watershed characteristics, which are the observable or the in- 
ferable characteristics determined from remote sensing, were used 
with the optimum model parameters to perform a correlation analysis 
(Task 6) and thereby determine statistical relationships between optimum 
model parameters and observable physical characteristics of the watershed. 

Five watersheds were reserved and used to test the relationships 
developed in the calibration and correlation tasks . Observable charac- 
teristics were used, together with the relationships determined by Task 
6, to generate sets of simulation model parameters. For each of these 
watersheds, the simulation model was then run, using actual climatolog- 
ical data as input. The simulated streamflow was then compared to the 
actual historical streamflow, to yield an evaluation or validation (Task 7) 
of the relationships. 


3-1 



3-2 



Figure 13. Feasibility Study Flow 




Task 8, Resolution/Ac curacy Analysis, was to be a determination 
of the performance characteristics of remote sensors needed to make 
application of the techniques practicable. The study was not able to ad- 
dress this task within the originally planned budget and schedule, but a 
very similar task is to be performed under another contract. 

This report is the result of Task 9 . 
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3 . 1 PHYSIOGRAPHIC DATA BASE 


Operation of the mathematical models (both OPSET and KWM) used 
in the hydrological simulation activity is governed by certain input para- 
meters which are related to the physical nature of the watersheds being 
simulated. These parameters must be measured or in some instances 
inferred from existing documentation such as topographical maps, aerial 
photographs, soil catalogs, survey reports, and other a priori information. 
This collected documentation, together with certain measured or inferred 
information therefrom, constitutes the physiographical data base whose 
construction is the objective of Task 1. Figure 14 depicts the activities 
involved in establishing the physiographic data base. Table 2 is a summary 
of the applicable documentation presently available in IBM files, and Table 
3 is a physiographic data base status summary for the 15 watersheds utilized 
in the feasibility study. 

The physiographic data required by the mathematical models is a 
series of parameters as listed and defined in Table 4. The methodology 
employed to derive these parameters from the physiographic documenta- 
tion is discussed in the following paragraphs. 

3.1.1 AREA 

The first step in determining the area of any given watershed is to 
locate on a topographical map or photomosaic the streamflow gage asso- 
ciated with that particular watershed. Then beginning at the stream gage 
location, an outline of the watershed is manually drawn by following con- 
tour lines and ridges until the area perimeter closes on itself back at 
the stream gage. Visualization of ridges, slopes, and drainage patterns 
is considerably enhanced by accentuating the waterways associated with 
a particular watershed. Accentuation is easily accomplished by darken- 
ing the rivers and streams on the map with a dark blue pencil. Once the 
watershed perimeter is determined, the area can be measured with a 
planimeter. 

3.1.2 OVERLAND FLOW SURFACE LENGTH (OFSL) 

The overland flow surface length is defined as the average distance 
excess precipitation, or runoff, must travel from ground impact until it 
reaches a permanent water course such as a river, stream, creek, branch, 
or crevasse. The average value of OFSL is determined by selecting 10 to 
20 points at random throughout the watershed as outlined on the topograph- 
ical map. A line perpendicular to the contour lines is then drawn from the 
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Figure 14. Physiographic Data Base 






Table 2. "Raw" Physiographic Data Acquired 
(Tennessee Valley Region) 


146 TOPOGRAPHIC MAPS, 7.5' QUADRANGLE 

385 AERIAL PHOTOGRAPHS, USDA B&W; 1:15,000, 1:24,000, and 1:36,000 

45 AERIAL PHOTOGRAPHS, NASA RB-57F, B&W PRINTS FROM IR NEGATIVES, 

1:120,000 

250 GEOLOGIC AND MINERAL RESOURCE SUMMARIES 


MISCELLANEOUS SPECIAL REPORTS AND DOCUMENTS 
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Table 3. Physiographic Data Base Status 


WATERSHED 

MAP 

COVERAGE 

PHOTO 

COVERAGE 

PARAMETERS 

QUANTIFIED* 

REMARKS 

WHITE HOLLOW 



NONE 

12/12 

PHOTOS REQUIRED TO 

LITTLE CHESTUEE 



PARTIAL 

12/12 

ESTIMATE VALUES FOR: 

S. FORK MILLS RIVER 



PARTIAL 

7/12 

(1) FIMP 

NOLAND CREEK 



PARTIAL 

7/12 

(2) VINTMR 

W. FORK PIGEON RIVER 



PARTIAL 

7/12 

(3) OFMN 

BIG BIGBY 



PARTIAL 

7/12 

(4) OFMNIS 

BIG ROCK 

100% 

PARTIAL 

7/12 


LITTLE RIVER 



PARTIAL 

7/12 

ALL OTHER PHYSICAL 

TRACE CREEK 



NONE 

7/12 

PARAMETERS (EXCEPT 

PINEY CREEK 



NONE 

7/12 

CHANNEL CAPACITY— 

RUTHERFORD CREEK 



NONE 

7/12 

CHCAP) CAN BE 

EMORY RIVER 



NONE 

12/12 

DETERMINED OR 

SEWEE CREEK 



NONE 

7/12 

INFERRED FROM 

TOWN CREEK 

1 

f 

NONE 

7/12 

TOPOGRAPHIC MAPS. 


Expressed as a fraction of the total number of parameters to be quantified. 
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Table 4. 


Physiographic Data Parameters 


AREA 

OFSL 

OFSS 

TAH 

FIMP 

FWTR 

VINTMR 

OFMN 

OFMNIS 

GWETF 

SUBWF 


Area of the Watershed in Square Miles 
Average Overland Flow Surface Length in Feet 
Average Overland Flow Surface Slope 
Time Area Histogram Parameters 

Fraction of Total Watershed Area Covered by Impervious Surface 
Fraction of Total Watershed Area Covered by Water Surface 
Vegetative Interception Maximum Rate in Inches/Hour 
Overland Flow Manning’s “N” for Pervious Surfaces 
Overland Flow Manning’s “N” for Impervious Surfaces 
Ground Water Evapotranspiration Factor 
Subsurface Water Factor 


(The significance and methods of quantification of these parameters are described in the text.) 



selected point to the nearest water course. The length of that line is then 
measured and tabulated. This same procedure is repeated for all selected 
points and the average distance of the 10 to 20 points is calculated and re- 
corded as OFSL. 

3.1.3 OVERLAND FLOW SURFACE SLOPE (OFSS) 

The OFSS is the average slope of the land contained within the 
boundaries of the watershed and is determined as follows. For con- 
venience, the same points selected for OFSL may be used to determine 
OFSS. From each point selected, the height differential from that point 
to the nearest watercourse can be determined by counting the contour 
lines and multiplying by the contour scale. Since the length of the over- 
land flow surface at that point was measured when determining OFSL, 

OFSS is then equal to the change in height (Ah) divided by the average 
of the local OFSS values . 

Table 5 is a sample worksheet showing the final average OFSS and 
OFSL values for a typical watershed. 

3. 1. 4 TIME AREA HISTOGRAM 

A time area histogram requires the division of the watershed into 
"n" subareas where each area is defined as a function of drainage time. 
Boundaries of subareas are known as isochrones, which may represent 
15- minute or one-hour intervals depending on the size of the watershed. 

A minimum of four isochrones is desired for operation in the simulation 
model. The steps required to develop a time area histogram are as 
follows : 

1. Utilizing the topographical map of the watershed, 
determine the longest watercourse in units of feet 
from the stream gage location to the furthermost 
point of the defined watershed. This will be noted 
as a stream length (L). 
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Table 5. Worksheet for OFSS/OFSL (Typical) 


Site 

Overland 

Altitude 

Local 

Selected 

Distance 

Change 

Slope 

(No.) 

(Ft.) 

(Ft.) 

(AH/S) 

1 

1,000 

160 

0.160 

2 

900 

100 

0.110 

3 

1,400 

260 

0.186 

4 

1,200 

120 

0.100 

5 

2,000 

120 

0.060 

6 

1,100 

80 

0.073 

7 

1,500 

160 

0.107 

8 

700 ■ 

160 

0.229 

9 

2,000 

180 

0.090 

10 

700 

180 

0.258 


12,500 


1.333 


OFSL = 12,500 -f 10 = 1,250 ft. 
OFSS = 1.333 -M0 55 0.1333 


NOTE: Above values are actuals for Rutherford Creek near Carters Creek, Tennessee 



2. Determine the change in elevation in feet between the stream 
gage location and the furthermost point in the basin as 
established in Step 1 . This change will be noted as Ah. 

3. Calculate the concentration time (T c ) in minutes utilizing 

the empirical equation 1 ^ 



where L is stream length and Ah is altitude difference 
‘ as defined previously. Concentration time is the time it 
takes for a drop of water to reach the gaging station from 
the watershed extreme perimeter. If L and Ah are measured 
in meters rather than feet, then the coefficient should be 
0.0195 rather than 0.0078. When T c is equal to or greater 
than 240 minutes, the isochrone intervals should be one hour. 
This will provide the minimum subareas required to operate 
in the simulation models. When T c is less than 240 minutes, 
then the isochrone interval should be 15 minutes. Since it is 
very unlikely that T c will calculate in exact multiples of 15 
minutes, the concentration time should be arbitrarily adjusted 
' to the nearest number divisible by 15. 

4. Calculate the number of subareas or isochrones by dividing 
the concentration time by either 15 minutes or 60 minutes. 

5. Calculate the average stream velocity (V) by dividing the 
stream length (L) by the concentration time (T c ). 



6. Calculate the average overland flow distance per isochrone by 
multiplying V by either 15 minutes or 60 minutes. 


3-11 



NOTE: Figure 15' shows an example of the calculations used to establish 
a time area histogram. 

Once the average overland flow distance per isochrone has been 
determined, the actual subdivision of the subdivision of the watershed 
can be accomplished as follows: 

1 . Beginning at the stream gage location on the topographical 
map and using a map measuring meter, track several water- 
courses from the gage station inland to a distance equal to 
the average overland flow length per isochrone for that 
particular watershed. Measured distances must always 
follow the watercourses. Several measurements should be 
made radially from the stream gage and a mark placed at 
each overland flow distance milestone. A line connecting 
these marks should then be drawn from perimeter to 
perimeter and this area identified as Isochrone 1. This 
isochrone identifies that area of the basin which would drain 
within the first 15 minutes (or one hour) after a storm. Then 
measuring from the first isochrone line, repeat the process 
to establish the second isochrone area. This is the area that 
would be drained from 15 to 30 minutes after a rain. Continue 
the process until the entire basin is subdivided into "n" areas. 
Figure 16 is an illustration of a typical watershed so divided. 

2. Once the watershed is subdivided, measure with a planimeter 
the area of each isochrone. These areas should be tabulated 
as both absolute magnitude in square miles and also as a 
fractional portion of the entire watershed. Table 6 is a typical 
data sheet with this information. Each isochrone area and 
fractional portion value is converted to punch cards and 
entered into the simulation model. 

3.1.5 IMPERVIOUS SURFACE FRACTION (FIMP) 

This is a fraction of the total watershed area covered by pavings, 
rooftops, streets, etc., whose runoff contributes directly to a stream. 
This value is usually estimated from aerial photography. The value is 
usually zero for rural areas except where there are large areas of 
exposed rock adjacent to watercourses . 
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Maximum Stream Length (measured): 
Maximum Attitude Drop — H(Map Scale): 


110,000 Ft. 
480 Ft. 



= 482 Minutes = 480 Minutes 

^isochrones = 480 + 60 = 8 

V = 110,000 Ft./480 Minutes = 229 Ft./Minutes 

Average Overland Flow Distance/Isochrone: 229 Ft./Min. x 60 Minutes = 13,750 Feet 

NOTE: Above calculations are actuals for Sewee Creek Watershed near Decatur, Tennessee. 


Figure 15. Time Area Histogram Calculation (Typical) 



Area: 

Concentration Time: 
Average Velocity: 
Maximum Stream Length: 
Delta Elevation: 


35.7 km 2 
75 Min. 

180 m/Min. 
13,500 Meters 
1,325 Meters 


1 


2 


3 


4 


WATERSHED 
BOUNDARY ' 



C. ) Divided by Isochrones 
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Table 6. Time Area Histogram Worksheet (Typical) 


Isochrone No. 

Area (Sq. Miles) 

Fraction of Total 

1 

3.9 

0.034 

2 

16.9 

0.144 

3 

6.8 

0.058 

4 

11.0 

0.094 

5 

18.3 

0.156 

6 

23.2 

0.198 

7 

15.3 

0.131 

8 

21.6 

0.185 


117.0 

1.000 


NOTE: Above values are actuals for Sewee Creek Watershed near Decatur, Tennessee 



3.1.6 WATER SURFACE FRACTION (FWTR) 


This is a fraction of the total watershed area covered by water 
surfaces at normal low flow conditions. The value is estimated from 
aerial photos and is virtually zero for watersheds containing neither 
lakes nor swamps. 

3. 1.7 VEGETATIVE INTERCEPTION MAXIMUM RATE (VINT MR) 

VINTMR is the maximum rate of rainfall interception by the water- 
shed vegetation expressed in inches per hour. The values range from 
0. 10 for grasslands to 0. 20 for heavy forest and is estimated from land 
cover interpretation from aerial photographs. A weighted average is 
used for watersheds having more than one type of vegetative cover. 

See reference 13. 

3.1.8 OVERLAND FLOW MANNING'S M N" FOR PERVIOUS SURFACES 
(OFMN) 

This is a roughness coefficient for overland flow derived from 
published tables dependent on estimated vegetative cover and soil 
usage. Weighted averages are used where different types of cover are 
in evidence. Values range from 0.018 for smooth earth to 0.100 for 
heavy forest. (See Table 7 for more details.) 

3.1.9 OVERLAND FLOW MANNING'S "N" FOR IMPERVIOUS SURFACES 
(OFMNIS) 

Derivation is the same as for OFMN but applies to impervious surfaces. 
Values are derived from published tables and vary from 0.013 for smooth 
asphalt to 0.017 for unfinished concrete. (See Table 7.) 

3.1.10 GROUND WATER EVAPOTRANSPIRATION FACTOR (GWETF) 

GWETF is a factor used to estimate the current rate at which 
swamp vegetation draws water from below the water table. This value 
is usually zero unless aerial photographs reveal a significant area of 
swamps within the watershed perimeter. In this case an appropriate 
value may be estimated using methods analogous to that used in 
determining FIMP. 


3-16 



3-17 


Table 7. 


Manning’s Roughness Coefficient for Overland 
Flow for Various Surface Types 


Watershed Surface 

Manning's "N" 

Smooth Asphalt 

0.013 

Concrete (Trowel Finish} 

0.013 

Rough Asphalt 

0.016 

Concrete (Unfinished) 

0.017 

Smooth Earth 

0.018 

Firm Gravel 

0.020 

Cemented Rubble Masonry 

0.025 

Pasture (Short Grass) 

0.030 

Pasture (High Grass) 

0.035 

Cultivated Area (Row Crops) 

0.035 

Cultivated Area (Field Crops) 

0.040 

Scattered Brush, Heavy Weeds 

0.050 

Light Brush and Trees (Winter) 

0.050 

Light Brush and Trees (Summer) 

0.060 

Dense Brush (Winter) 

0.070 

Dense Brush (Summer) 

0.100 

Heavy Timber 

0.100 


SOURCE: Reference 15, pp. 110— 113. 




3.1.11 SUBSURFACE WATER FACTOR (SUBWF) 


SUBWF is the fraction of moisture entering groundwater storage 
which leaves the basin through subsurface flow (underground rivers, 
subterranean caverns) not measured by the stream gage. This value is 
usually zero unless such subterranean geophysical anomalies are known 
to exist. 


3.2 HISTORICAL DATA BASE 


The historical data required by the simulation and calibration 
models is similar in format to the physiographic data described in 
paragraph 3.1. The basic difference is that while all physiographic 
data is either measured or inferred from observations, historical data 
is converted from existing written records to a digital format for input 
to the models. Where precipitation gages are not uniformly distributed 
about a given watershed and are a combination of hourly and daily gages, 
Thiessen analysis is used to derive "weighting factors" to apply to 
the historical data. 


The historical data base generated by IBM for this hydrology 
application study is constructed primarily from five different specific 


types of 

data: 

1 . 

Precipitation data 

2. 

Stream stage charts 

3. 

Stream gage rating tables 

4. 

Daily discharge data 

5. 

Evaporation data 


The conversion of this data to the digital forms usable by the simulation 
and calibration models is represented in Figure 17. Table 8 is a listing 
of the historical data collected and presently on file in IBM, and Table 9 
is a summary of the utilization of that data base as applied to the 
initial 15 watersheds being studied. The following paragraphs describe 
in some detail the procedures used to extract the desired information from 
the historical data base. 
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Table 8. "Raw" Historical Data Acquired 
(Tennessee Valley Region) 

NWS DAILY CLIMATOLOGICAL DATA TABLES, 1948-70 

TVA DAILY PRECIPITATION TABLES, 1948-70 

NWS HOURLY PRECIPITATION DATA TAPES (5 REELS), 1960-70 

NWS DAILY PRECIPITATION DATA TAPES (9 REELS), 1963-70 

TVA DAILY PRECIPITATION DATA TAPE (1 REEL), 1968-71 

USGS DAILY SURFACE WATER RECORDS, 1961-70 

TVA DAILY DISCHARGE TABLES, 1940-70 

USGS DAILY DISCHARGE DATA TAPE (1 REEL), 1939-70 


SELECTED STREAM STAGE RECORDER CHARTS 
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Table 9. Processed Historical Data 


MEAN BASIN/HOURLY PRECIPITATION 

• 15 MERGED RECORDS 

• 100 WATER YEARS 

• 294 WATER YEARS INPUT DATA 


DAILY DISCHARGE 

• 15 RECORDS 

• 100 WATER YEARS 


STORM EVENT HYDROGRAPHS 

• PEAK DAY, HOUR, MAGNITUDE 

• 500 SIGNIFICANT STORMS 


EVAPORATION 

• 15 RECORDS 

• 100 WATER YEARS 



3. 2. 1 DATA AVAILABILITY SUMMARY 


Precipitation data is gathered throughout the Tennessee Valley by 
many precipitation stations which are owned and/or managed by either 
TVA, the National Weather Service (NWS), or private corporations. In 
all cases the data is recorded either hourly or daily and the information 
stored on magnetic tape and/or written documents. Magnetic tape 
storage is limited primarily to data collected from certain stations since 
1960; therefore, most precipitation data of value to the study is recorded 
in tables. The following paragraphs detail the step-by-step procedures 
necessary to obtain the hourly precipitation data required by the models 
for a given watershed. 

Once a subject watershed has been identified, several initial 
questions must be answered: 

1 . How many precipitation stations exist within the boundaries 
of the watershed? 

2. How many precipitation stations are closely adjacent to 
the watershed and how close are they? 

3. How many of the stations are hourly and how many are daily? 

4. What are the years of operation for the stations? 

5. How much of the existing historical data exists as written 
tables and how much on magnetic tape? 

The answers to the above questions exist in different reports, lists, maps, 
documents, etc. , and require considerable manual effort to compile. 

Figure 18 is a sample of a matrix type data worksheet developed to 
correlate the answers to the questions. The following paragraphs explain 
the notation used in Figure 18. 

The watershed name, approximate latitude and longitude boundaries, 
centroid, and station names are self-explanatory. The station descriptions 
are (1) S.G. for stream gage (required to determine the operational years 
to be selected for input to the simulation model), (2) PG-H and PG-D for 
precipitation gage-hourly and daily, respectively, and (3) E.S. for 
evaporation station. 
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WATERSHED NAME- NOLAND CREEK NEAR BRYSON CITY 

NORTH CAROLINA (13.8 SQ. MILES) 



WATERSHED BOUNDARIES (APPROX.) LATITUDE 3528 

TO 

3534 

LONGITUDE 8327 

TO 

8331 



WATERSHED CENTROID (APPROX.) LATITUDE __ 


3531 


.... lONRlTlinF 

8329 

— 


STATION NAME 

Station 

TVA 

NWS 

LOCATION 

Read 

— 

Data 

OPERATIONAL YEARS 

De&cr. 

No. 

No. 

Lat. 

Long. . 

Time 

Type 

49 

50 

51 

52 

i 53 

54 

155 

U 

57 

58 

59 


61 


63 

\u 


66 



69 

U 

NOLAND CHEEK 

S.G 

(3-5135) 
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8330 

R 
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Hourly Discharge Data Available. 
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Figure L8 . Historical Data Availability 



The station latitude and longitude are important in determining the 
proximity of the precipitation gages to the watershed. The location of the 
various gages is required for the merging and synchronizing process 
described in paragraph 3.2.2. The read time is either "R M for continuous 
recording or is specified for daily stations (e.g. , 2400 means the gage is 
read daily at midnight). 

The data type is either "T" meaning magnetic tape and/or "D M for 
written documentation. In those instances where both types are indicated, 
the upper line of the adjacent bar chart represents those operational years 
where the data exist on magnetic tape. The lower line of the bar chart 
reflects the years that the station has been in operation and for which 
precipitation data (tape and/or table) is available. The ’’hourly discharge 
data available" entry at the bottom of the charts represents those years 
for which strip chart hydrographs are available. 

Once this worksheet is completed, the water years to be used for 
that particular basin can be selected, utilizing available stream gage 
data and the best of the available precipitation gages . After selection 
of operational water years and precipitation stations, conversion of the 
raw data to the digital data required by the models can begin. 

3.2.2 HOURLY PRECIPITATION DATA 

Hourly precipitation data in digital form is the primary input to 
KWM and OPSET. In a very small watershed having its own hourly 
precipitation gage (e.g. , White Hollow) one can with reasonable safety 
assume that the gage reading applies uniformly to the entire watershed. 
This assumption (which is implicit in both programs) departs from reality 
more and more with increase in watershed size. It has been necessary, 
as part of the study, to implement a method whereby several precipitation 
records are used to synthesize a single hourly rainfall history for each 
basin. 


The number of precipitation stations associated with any given 
watershed may vary from one station located 20 or 30 miles from the 
watershed centroid to 5 or 6 stations located within or closely adjacent 
to the watershed boundaries. 


3-24 



Typically, a watershed will have one or two hourly stations, and 
one or more daily stations. In addition to the varying distances of these 
stations from the centroid, the reading time for the daily stations might 
be different. It is also quite likely that data will appear from the several 
gages in both magnetic tape and tabular formats. The latter must be 
manually extracted from the tables and converted to punched data card 
format. 

The precipitation gage outputs are assigned weighting factors, using 
the Thiessen technique 2 ? 15, in accordance with their physical locations 
relative to the basin centroid. A software program, designated MERGE 
AND SYNC in Figure 17 which was developed by IBM Huntsville, auto- 
matically performs the interpolation and correlation of the precipitation 
data. This program accepts all precipitation data, the reading time for 
each daily station, and the weighting factor developed from the Thiessen 
Analysis, and produces an hourly precipitation record for the applicable 
water years associated with a given watershed. This hourly precipitation 
data record is then used as one of the climatological inputs required by 
the models. 

3.2.3 SELECTED STORM DATA 

For operation of the OPSET program it is necessary to select up to 
five flood hydrographs for each of the years for which the model is to be 
calibrated. This requires a manual search of precipitation and discharge 
records to select storms useful to the calibration. The digitized input 
data include the number of hydrographs chosen and three parameters 
related to each hydrograph: day of occurrence of the flood peak, hour of 
occurrence of the flood peak, and flow rate at the peak. These hydrographs 
parameters are essential for the OPSET program to determine watershed 
model routing parameters, so that total flows will represent accurate 
predictions, with respect to the time of occurrence of hydrograph peaks 
as well as the total volume of flow for a given period of time. In practice 
the selected storm hydrograph parameters are not available in daily 
discharge records. It is necessary to obtain them from the strip charts 
produced by the stream gage recorders. Rating tables are also digitized 
and stored for conversion of gage height readings into flow rate. 

The procedure employed to obtain this data requires manual analysis 
of each strip chart and manual recording of the rise and fall of the stream 
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gage on an hourly basis. The time frame should extend from midnight 
of the day in which the storm occurred until some time at which the 
stream height returns to or approaches its initial stage. This hourly 
height recording is then formatted for entry into the computer where a 
subroutine will fetch the appropriate rating table into memory and 
convert the data to cubic feet per second. This flood hydrograph data 
is then in a usable form when required by the simulation model. 

3.2.4 DAILY DISCHARGE DATA 

Daily discharge data is the volume in cubic feet of water per day 
that flows past any given stream gage. This data exists on magnetic 
tape and/or written tables for all stream gages in the Tennessee Valley. 

The data format which exists on magnetic tape must be altered to be 
compatible with the simulation model. Where the data exists in written 
tables, it is necessary to manually extract that information, convert to 
punched card format, and develop a listing compatible with model requirements. 

3.2.5 EVAPORATION TABLES 

Evaporation data appear in Climatological Data publications of the 
National Weather Service. Unfortunately, the number of pan evaporation 
stations is too limited to provide complete coverage. The nearest evapo- 
ration station may be as much as 100 miles from the watershed. Additionally, 
the station may be associated with a large lake or reservoir which has 
evaporation rates different from those of an interior watershed in a pre- 
dominantly mountainous region. Preparation of the evaporation data is 
similar to that for daily discharge data in that the rates and pan evapo- 
ration coefficients are read from published tables, punched onto cards, 

and a computer- compatible listing generated for the identified watershed. 

* 

3.3 MASTER WATERSHED DATA BANK 


The digital products of Task 1 (Physiographic Data Base) and Task 2 
(Historical Data Base) are stored in a Master Watershed Data Bank, which 
is in a digital tape format, along with control options and other logical 
inputs needed for the operation of OPSET and KWM as indicated in Figure 19 . 
This tape-based data bank provides flexibility in operating the models and 
obviates storage of some 1,200,000 punched cards. 
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Figure 19. Task 3: Master Watershed Data Bank 



Considerable effort has been devoted to acquiring and formatting 
the raw data and creating this master data bank. It is a by-product of 
the study which is potentially valuable to others engaged in hydrologic 
or related research using the same geographic area. 

3.4 CALIBRATION (TASK 4) 

3.4.1 USEOFOPSET 

The OPSET model is used to calibrate a set of watersheds chosen 
from a region, generally using three water years of historical data for 
each basin, determining a set of model parameters from the calibration 
runs, and testing the accuracy of the calibration using the KWM. Cali- 
bration of several watersheds using OPSET generated sets of model 
parameters that produced accurate simulations on a daily flow and 
monthly flow basis which unfortunately did not accurately simulate the 
low flows. The subroutine that estimates the Base Flow Recession 
Constant (BFRC) has been modified, and results indicate an improvement 
in accuracy of simulation of low flows. Some manual "fine tuning" has 
been used to obtain a set of model parameters that represents the water- 
shed very well. 

Figure 21 shows the activities involved in Task 4 (Calibration). Data 
pertaining to the watershed to be calibrated are fed into the OPSET program, 
which is then run to estimate a set of model parameters. This step gets the 
task "into the ballpark. " A simulation is then run using KWM and the IBM 
analysis/evaluation routines. Simulation accuracy is evaluated with 
respect to total annual runoff, monthly flow, daily flow, statistical indices, 
and selected storm hydrograph characteristics. Based on these evaluations, 
parameters are adjusted, and a new simulation run, followed by another 
evaluation. This process goes through several iterations until simulated 
flow matches observed flow with acceptable accuracy in all criteria of 
interest to the analyst. The choice of parameters to adjust, direction 
and magnitude of the adjustments depend upon the judgment of the analyst. 
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Figure 20. Task 4: Calibration 




Sensitivity analyses performed to date have produced invaluable guidance 
to the manual- adjustment activity, reducing the subjectivity and eliminating 
the requirement that the analyst be skilled in hydrology. 

3. 4. 2 PARAMETERS ESTIMATED BY CALIBRATION 

The calibration process (a combination of automatic and manual 
adjustments as summarized above) was used to quantify the following 
simulation model parameters for each of ten watersheds. See Reference 
13 for more complete description. 

1. BFRC, Base Flow Recession Constant, governs the rate 
at which groundwater flow recedes in the model. 

2. IFRC, Interflow Recession Constant, governs interflow 
recession. 

3. BUZC, Basic Upper Zone Capacity, is an index for 
estimating the capacity of the soil surface (upper zone) 
to store water in interception and depression storage. 

4. SUZC, Seasonal Upper Zone capacity adjustment constant, 
is used to adjust upper zone variations in vegetation and 
cultivation. 

5. LZC, Lower Zone Capacity, is an estimate of the capacity 
of the basin soil to hold water. Decreasing LZC in the 
model has the effect of increasing synthesized runoff. 

6. BMIR, Basic Maximum Infiltration Rate, is the index used 
to control the basic rate of moisture inflitration. This is 

a parameter to which simulation accuracy is very sensitive, 
particularly as it affects storm peaks. 

7. SIAC, the Seasonal Infiltration Adjustment Constant, is 

an evaporation-infiltration factor relating infiltration rates 
to evaporation rates to account for more rapid infiltration 
during warmer periods. 

8. ETLF, Evapotranspiration Loss Factor, is an index used 
to estimate the maximum rate of evapotranspiration which 
could accur within the watershed under current conditions 
of soil moisture content. 

9. BIVF, Basic Interflow Volume Factor, controls time distri- 
bution and quantities of moisture entering interflow. Increasing 
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BIVF tends to reduce storm runoff peaks and extend 
hydrograph recession limbs. 

10. NCTRI, Number of Current Time Routing Increments, is 
the number of subareas into which the basin should be 
divided, given 15 minute or one hour separation between 
isochrones. 

11. CSRX, Channel Storage Routing Index, is used to account 
for channel storage when channel flows are less than half 
of channel capacity (CHCAP). Channel storage effects 
are simulated by having the hydrograph time routed to the 
mouth of the watershed through an imaginary reservoir. 

12. FSRX, Flood plain Storage Routing Index, is used to account 
for channel storage plus flood plain storage when stream- 
flows are greater than twice the channel capacity. Between 
one-half and twice channel capacity, the program inter- 
polate values between FSRX and CSRX . 

13. CHCAP, Channel Capacity, is that value of streamflow, 
measured at the gage, at which a transistor is made from 
channel routing to flood plain routing. In mountainous 
watersheds, this is not an oritical parameter, and OPSET 
seldom adjusts it. 

After OPSET has adjusted the parameters listed above to achieve a 
”best match” based on mean daily streamflow, it is usually found that 
synthesized flood peaks fail to match observed peaks, in magnitude and/or 
time. Since the study attempted to address as wide a variety of applications 
as practicable, some manual ’’fine tuning” was undertaken to achieve 
an acceptable match between synthesized and observed flood peaks while 
maintaining an acceptable correlation between synthesized and observed 
mean daily and monthly flows. This manual adjustment process requires 
some knowledge of the hydrologic processes occurring in the watershed 
and some subjective judgment. No firm rules or recipes have been 
developed, but the following are useful guidelines. 

1. Overall results are affected by soil moisture capacities 
and infiltration rates (LZC, BUZC, BMIR) and their related 
seasonal adjustment constants (SUZC, SIAC). 

2. Initial storages (LZS and UZS) can be varied to improve 
accuracy in the first two months of the multiyear simula- 
tion. 
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3. Summer storm peaks are affected more than winter 
storm peaks by changes in SIAC and SUZC, and the 
latter has more influence on mean daily flow in drier 
months than the former. 

4. Consistent phasing errors (differences in times of occurrence 
between synthesized and observed flood peaks) can usually 

be reduced by adjustment if the number and sizes of 
subareas in the time-area histogram. Phasing errors 
which appear random are attributable to errors in input 
precipitation and/or evaporation data. 

5. Since there are parameter interactions, all performance 
indices should be re- checked after any parameter adjust- 
ment and others re- adjusted as needed until a "best simula- 
tion" is realized. 

3.4.3 CALIBRATION RESULTS 

The calibration process was applied to ten Tennessee Valley 
watersheds to derive the sets of parameters shown in Table 10. Parameter 
definitions are listed in Table 11. Each set is the one which produces 
simulations that are acceptable with respect to all indices of performance, 
as shown in Table 12. The mean monthly and daily streamflows show 
good correlation coefficients. 

From each multi-year simulation, the largest winter storm was 
included in the table, even though others showing more felicitous results 
could have been chosen. For some of the large winter storms (e.g., 

Sewee Creek, 3-13-63), simulated flood peaks vary from the observed 
peaks by factors in magnitude and several hours in time. Experience 
in the study shows that (1) input data are often erroneous and (2) good 
simulation fidelity with respect to storm/flood events should not be 
expected in a free- running or "open loop" model. In an operational 
application, moisture storage conditions and simulated flow rate would 
be periodically adjusted in accordance with actual observed conditions and 
in a short-term forecast (24 to 72 hours), a prediction accuracy within 
five per cent could confidently be expected. 

Several plots were generated early in the study after the first cali- 
brations were completed. They were useful at first, but more time- 
consuming than the information they yielded would justify. It was found that 
simulated and observed storm events could be compared effectively, for 
evaluation purposes, on the basis of peak flow magnitude and time. Some 
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Table 10. Simulation Model Parameters After Calibration 


WATERSHED 

WATERSHED 

PARAMETERS 

PARAMETERS QUANTIFIED BY CALIBRATION 

SOIL 

CHARAC- 

TERISTICS 

SOIL MOISTURE | 

CHANNEL ROUTING 

AND GROUNDWATER 



LAND F 

LOW 


NO. 

NAME 

AREA 



FIMP 

FWTR 

BMIR 

LZC 

ETLF 

SUZC 

S1AC 

BUZC 

BIVF 

CSRX 

FSRX 

BFRC 

BFNLR 

OFSL 


OFMN 

IFRC 

OFMNIS 

AWC 

PERM-A 

1 

WHITE HOLLOW 

2.68 

.001 

.001 

9.0 

6.0 

.20 

.50 

.40 

.90 

.50 

.99 

.97 

.98 

.85 

1000 

.30 

.10 

.10 

.015 

3.71 

1.70 

2 

LITTLE CHESTUEE 

8.26 

.013 

0.0 

10.0 

6.0 

.20 

.20 

.40 

.50 

0.0 

.93 

.97 

.99 

.85 

630 

.237 

.061 

.10 

.002 

8.75 

1.20 

3 

SOUTH FORK MILLS 

9.99 

.001 

0.0 

12.9 

7.6 

.07 

.20 

.20 

.20 

.58 

.96 

.97 

.99 

.85 

1140 

.315 

.05 

.10 

.015 

7.90 

4.54 

5 

WEST FORK PIGEON 

27.6 

.001 

.008 

4.0 

7.0 

.10 

.15 

.15 

.15 

0.0 

.95 

.95 

.99 

.85 

1380 

.413 

.05 


.015 

8.82 

4.0 

7 

BIG ROCK CREEK 

24.9 

.033 

.005 

5.0 

4.0 

.07 

.15 

.15 

.40 

0.0 

.94 

.95 

.70 

.50 

980 

.155 

.04 

.10 

.015 

5.30 

.77 

9 

TRACE CREEK 

30.4 

.039 

.005 

6.5 

4.6 

.11 

.70 

.90 

.39 

0.0 

.97 

.97 

.96 

.85 

850 

.129 

.05 

.10 

.015 

1.80 

1.43 

12 

EMORY CREEK 

j 

83.2 

.005 

.001 

7.0 

3.9 

.10 

.20 

.20 

.30 

.78 

.94 

.94 

.84 

.85 

1270 

.298 

.07 

.10 

.010 

5.55 

2.67 

13 

1 

SEWEE CREEK 

117 

.001 

.001 

4.0 

3.5 

.10 

.25 

.15 

.25 

.50 

.98 

.99 

.99 

.85 

1060 

.165 

.05 

.10 

.015 

4.92 

.50 

14 

TOWN CREEK 

141 

.002 

.001 

7.0 

4.0 

.25 

.25 

.30 

.25 

.50 

.94 

.90 

.93 

.85 

1550 

.062 

.05 

.10 

.014 

6.19 

4.3 

15 

POPLAR CREEK 

82.5 

.010 

.001 

4.0 

3.0 

.15 

.20 

.15 

.15 

1.0 

.97 

.98 

.98 

.85 

1320 

.329 

.05 

.10 

.015 

2.02 

1.08 


NOTE: Watershed parameters and soil characteristics are included in the table to make a complete base for correlation. 
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Table 11. Model Parameter Definitions 


SYMBOL 

DEFINITION/SIGNIFICANCE 

BIVF 

Basic Interflow Volume Factor; Indexes distribution and quantities of moisture entering interflow. 

BFNLR 

Base Flow Nonlinear Recession Adjustment Factor 

BFRC 

Base Flow Recession Constant 

BM1R 

Basic Maximum Infiltration Rate; relates to "A" horizon permeability. 

BUZC 

Basic Upper Zone Storage Capacity 

CHCAP 

Channel Capacity Indexed to Basin Outlet; flow rate indicated at transition from channel 
routing to flood plain routing. 

CSRX 

Channel Storage Routing Index; used to account for channel storage when channel flows are less 
than half of CHCAP. 

ETLF 

Evapotranspiration Loss Factor; indexes moisture depletion rate through evapotranspiration. 

FSRX 

Flood Plain Storage Routing Index; used to account for channel plus flood plain storage when 
streamflows are greater than twice CHCAP. The model synthesizes a routing index when 
streamflow is between Vs^CHCAP and 2*CHCAP. 

IFRC 

Interflow Recession Constant; indexes interflow recession rate, effective only in conjunction 
with a non-zero BIVF. 

LZC 

Lower Zone Capacity; indexes lower zone soil moisture storage capacity; relates to plant- 
available water capacity. 

NBTRI 

Number of Base Time Routing Increments; used by the model in routing by Clark's method. 

SIAC 

Seasonal Infiltration Adjustment Constant; relates to changes in vegetative cover with changes 
in season. 

suzc 

Seasonal Upper Zone Storage Capacity Factor; BUZC and SUZC together regulate the capacity 
of the upper zone to hold moisture. 

I VINTMR 

Vegetative Interception Maximum Rate. 
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Table 12. Simulation Accuracy After Calibration 



IVATCDCUCn 



YEARS 

OF 

RECORDS 

STREAM FLOW 

LARGEST WINTER STORM j 


VHA 1 tnontU 

CORRELATION 

COEFFICIENT 

DAILY, CFS 

DATE 

PEAK I 

NO. 

SHORT NAME 

AR 

EA 

MEAN 

. MAX 

FLOW, CFS 

TIME, HR 

Ml 2 

KM 2 

YRS 

FROM 

THRU 

M'LY 

D'LY 

CBS 

SIM 

OBS 

SIM 

OBS 

SIM 

OBS 

SIM 

1 

WHITE HOLLOW 

2.7 

6.9 

7 

57 

63 

.95 

.91 

3.79 

3.69 

147 

78 

3—12—63 

270 

98 

6 

4 

2 

LITTLE CHESTUEE 

8.2 

21.3 

7 

50 

56 

.92 

.87 

13.0 

13.9 

388 

447 

4-15-56 

1,030 

895 

18 

18 

3 

SOUTH FORK MILLS 

9.9 

25.7 

3 

66 

68 

.84 

.80 

32.2 

32.5 

713 

706 

2-13-66 

1,618 

1,273 

6 

8 

5 

WEST FORK PIGEON 

27.6 

37.3 

9 

56 

64 

.93 

.82 

78.1 

79.1 

1600 

2641 

3-5-64 

4,714 

3,368 

25 

27 

7 

BIG ROCK CREEK 

24.9 

64.5 

7 

55 

61 

.95 

,93 

40.3 

40.9 

3800 

3895 

3-7-61 

6,518 

5,923 

24 

26 

9 

TRACE CREEK 

30.4 

78.7 

6 

64 

69 

.93 

.88 

1 

39.8 

i 

37.5 

2310 

1547 

12—4—64 

3,524 

2,579 

8 

8 

12 

EMORY RIVER 

83.2 

215 

4 

53 

56 

.95 

.83 

142 

145 

5460 

3150 

3-22-55 

13,000 

8,925 

3 

1 

13 

SEWEE CREEK 

117 

303 

9 

60 

68 

.93 

.84 

192 

142 

6450 

5658 

3—13—63 

20,400 

8,561 

35 

46 

14 

TOWN CREEK 

141 

365 

6 

61 

66 

.95 

.90 

282 

216 

9260 

8384 

3-25-64 

10,699 

9,233 

48 

49 

15 

POPLAR CREEK 

82.5 

214 

7 

62 

68 

.96 

.90 

157 

142 

6200 

4068 

3-12-63 

6,389 

5,252 

38 

31 




of the plots generated early in ftie study are shown in Figures 21 (whose 
time scale is unavoidably too small for good interpretation) through 25. 
Except for some graphical enhancement for the sake of good reproduction, 
all the plots except that of Figure 23 were generated automatically. 
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STREAMFLOW, CUBIC FEET/SECOND 



80 120 160 200 240 280 320 

DAYS FROM OCTOBER THROUGH SEPTEMBER 

Figure 21. One Year Simulated and Observed Flow, White Hollow, 1961 
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STREAMFLOW, LITERS/SECOND 



STREAMFLOW, CUBIC FEET/ SECOND 



Figure 22. One Month Simulated and Observed Flow, White Hollow, March 1961 


STREAMFLOW, LITERS/SECOND 








Figure 23. Selected Storm, November 16-23, 1957, White Hollow 


PRECIPITATION, INCHES 
(MM) 



MEAN DAILY OBSERVED AND SIMULATED STREAMFLOW (cfsl 



MARCH (DAYS) 

OBSERVED MONTHLY STREAMFLOW = 1301 CFS SIMULATED MONTHLY STREAMFLOW = 1156 CFS 


Figure 24. One Month Simulated and Observed Flows, Little Chestuee, 
March 1951 
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HOURLY OBSERVED AND SIMULATED STREAMFLOW (CFS) 



♦LITTLE CHESTUEE, TENN. 8.24 SM. SN5103 03/28/51-03/29/51 STUDY LC08 




OBSERVED RUNOFF = 2.25 INCHES TIME {HOUR) SIMULATED RUNOFF = 2.81 INCHES 


Figure 25. Simulated and Observed Runoff, Major Winter Storm, 
Little Chestuee, March 28-29, 1951 
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MEAN DAILY OBSERVED AND SIMULATED STREAMFLOW (CFS> 



Figure 26. Simulated and Observed Flows, Emory River, February 1956 
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3.5 SENSITIVITY ANALYSIS (TASK 5) 

-The success of a-general hydrologic model is measured by its 
ability to simulate streamflow sequences that match observed records. 

The model utilizes time sequences of climatological data, information 
on the physical watershed, and a set of values for model parameters. 
These model parameters theoretically relate to watershed physical 
characteristics; operationally, those values are estimated by a sequence 
of trials and adjustments ending in an acceptable match. The quality of 
a given trial is determined by the closeness with which the observed and 
simulated flows agree during every simulation period. The manifest 
impracticality of making all these comparisons in evaluating a trial 
simulation requires selection of a small number of quality indices. 

One such index is the total annual flow; however, a large number of 
combinations of parameter values will give the same total annual flow. 
Therefore we need to differentiate in selecting among these combinations 
by adding other indices such as mean, standard deviations, root sum 
squares, daily correlation coefficients, low and high flows, and hydro- 
graph characteristics. 

If one is to adjust a trial set of model parameter values in order 
to improve the matching of the observed and simulated streamflows, he 
needs information on the effect a given change in a given parameter value 
will have on the simulated streamflows. The universe of such information 
is a multidimensional response surface of values for each index at points 
representing each set of parameter values. The response surface of 
interest is bounded by the reasonable range for each individual parameter. 
Comprehensive mapping is impractical. 

The information on the direction and the rate of change of the index 
values per unit change of each parameter at the point representing the set 
is important and very useful. Theoretically an infinite number of points 
could be analyzed, but a good one to start with is one accepted after trial 
and error modeling to achieve a "best" match. 

The advantages of the sensitivity analysis lie in its ability to: (1) 
guide manual adjustments or programmed adjustments in a computerized 
optimization routine , indicating how much to change a given mode param- 
eter to effect a given correction in a statistical index and/or hydrograph 
characteristic, and also the effect it will have on other indices; (2) show 
which parameters need to be estimated carefully and which ones require 
only rough approximation; (3) indicate the critical parameters on which 
effort should be concentrated to find a correlation with watershed physical 
characteristics that can be observed by remote sensing; and (4) provide 
guidance on what difference in runoff characteristics to expect from a 
given difference in physical characteristics of two watersheds. 
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Sensitivity analysis is time consuming and costly because of the 
large number of computer runs required. Also, a completed sensitivity 
analysis may be invalidated by a change in the model, but recently the 
models have been stabilized sufficiently to make the analysis worthwhile. 
Figure 27 shows the interfaces of Task 5 (Sensitivity Analysis) with other 
study tasks. 

Two sensitivity analyses were performed, for the purpose of 
guiding manual calibrations, on two small basins: White Hollow and 
Little Chestuee. The results of the latter are tabulated in Table 13. 

The column headings represent various mean daily streamflow statis- 
tics and hydrograph characteristics that were selected to assess the close- 
ness of simulated flow to recorded flow. No single statistical quality 
index can be used in reckoning the quality of a match. Different users 
need to match different statistical indices depending on their objective, 
interest, or application (e.g. , flood control, low flows, water supply). 

None of the parameters except NBTRI affect flood peak timing. Flood 
peak times are therefore not included in Table 13. Parameter values 
producing the "best" simulation appear in parentheses in the PARAM 
column. 

The data shown in Table 13 need further interpretation to be in- 
formative. A "Unit Sensitivity”, defined by the following division 

US = (percent change in simulation result)/(percent change 
in parameter) 

was calculated for each parameter with respect to the simulation result 
of interest, as a measure of parameter effectiveness. The results for 
the Little Chestuee basin are summarized in Table 14. 

Examination of the results suggests a general scheme of parameter 
change that should improve accuracy of low flow simulation, and the peak 
of the biggest storm. Increasing the BFRC value from 0. 98 to 0. 99 would 
improve the summer low flow simulation. A reduction in ETLF from 0. 165 
to 0. 05 would increase the annual flow volume, and counteract the reduction 
due to BFRC increase. Reducing the CSRX value from 0. 97 to 0. 96 should 
increase the peaks of all floods and therefore better the peak of the biggest 
storm. Results of experiments with these combinations were used to advan- 
tage in the calibration task. 
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Figure 27. 


Task 5: Sensitivity Analysis 


TO TASKS 8, 9 
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Table 13 


Sensitivity Analysis of Little Chestuee (Part 1 of 4) 


SENSITIVITY ANALYSIS cnn LITTLE CHESTUEE WATERSHED WATER YEAR _19&J DRAINAGE AREA 8 . 24 Sfia a i e M iles 


♦ 

DAILY STREAMFLOW STATISTICAL SUMMARY 




STORM ANALYSIS SUMMARY | 


STORM OF 
1 / 14/51 

#1 

STORM OF 

STORM OF 

1 STORM OF 

FLOW 

STATISTICS 

ERROR 

STATISTICS 

AVERAGE 

ERRORS 

LOW FLOW 
STATISTICS 

2 / 1/51 

#2 

3 < 29 / 51 

#3 

9/ 22 > 51 

#4 

MEAN 

StD 

nt’v. 

ROOT 

SUM 

SQUARE 

OAfLY 

CORA, 

coEir. 

LOW 

FLOWS 

(4.5) 

HIGH 

FLOWS 

03.1) 

MEAN 

OF 

1 Aug, 1 

MEAN 

OF 

1 1 

LOWEST 

PEAK 

8/0 

UNI 

PEAK 

n/o 

PEAK 

R/ O 

peak 

R / O 

1 1 N | 



OF 

1 Feb. » 

(CFS) 

(HOUR) 

(CFS) 

<Houm 

( in | 

(CFS) 

IHOVJ a I 

1 m ) 

(CFS) 

HOOP' 


OBSERVED 

CONFIGURATION 

14.01 

27.9 





5.2 5 

7.24 

11.0 

191 

21 

0.49 

798 

9 

1.61 

954 

2 

2.25 

230 

15 

.20 

BEST SIMULATED 
CONFIGURATION 

14 70 

33 R? 

194 

.963 

♦0.7 

+1.6 

4.50 

10.97 

9.7 

222 

19 

0.61 

950 

8 

1.86 

954 

4 

2.B3 

231 

16 

.23 

1 D 

par am 

VALUE 

WATER 

5HED PAR 

AMETE 

is 


















51 

AREA 

4.12 

7.38 

16.83 

271 

.963 

-2.5 

-21.7 

2.24 

5.46 

4.9 

111 

19 

.30 

475 

8 

.03 

477 

4 

1.41 

115 

16 

.14 

52 

18.24) 

12.36 

22.17 

50.41 

498 

.963 

+3 9 

+24.9, 


16.46 

14.6 

333 

19 

.91 

1424 

8 

2. B0 

1432 

4 

,.JL£4_ 

346 

16 

.43 

53 

FIMP 

0.00 

14.59 

33.33 

183 

.963 

+0.6 

+0.2 

4.36 

10.71 

9.9 

’219 

19 

.60 

944 

8 

1.B5 

9B0 

_4 

—2,80 

216 

— 1E-. 

.,2E. 

5? 

(0 013) 

0.052 

15.36 

34.74 

216 

.959 

+0.9 

+5.8 

4,90 

11.76 

9.4 

232 

19 

.64 

966 

e 

1.91 

969 

4 

2.90 

273 

16 

.35 

56 


0.10 

16.07 

36.17 

249 

.952 

♦ 1.1 

+ 11.0 

5.33 

12.75 

9.1 

244 

19 

.68 

987 

8 

1.98 

987 

4 

2.99 

326 

16 

43 

59 


0.20 

17.56 

39 48 

330 

.930 

+ 1.7 

+21.8 

6.42 

14.81 

8.1 

269 

19 

.75 

1030 

B 

2.10 

1023 

A 

3.19 

437 

16 

,55 

60 


0.50 

22.04 

51.20 

612 

.652 

+3.2 

+54.0 

9.54 

21.15 

5.5 

343 

19 

.99 

1158 

8 

2.48 

1134 

4 

3.77 

775 

16 

1.09 

63 

FWTR 

0.01 

14.67 

33.95 

19B 

.962 

+0.4 

+2.5 

4.09 

10.71 

9.6 

224 

19 

.61 

954 

8 

1.88 

958 

4 

2.84 

240 

16 

.30 

65 

10.0) 

0.10 

14.30 

3644 

246 

.952 

-1.1 

+10.1 

3.34 

B.95 

6.2 

244 

19 

.67 

989 

8 

1.97 

989 

4 . 

2.99 

327 

16 

.42 

6? 


0.50 

18.20 

49.27 

557 

.862 

+04 

+44.3 

7.24 

15.27 

2.4 

332 

19 

.93 

1144 

8 

2.41 

112B 

4 

3.67 

717 

16 

.98 




SOIL MO 

ISTURE P 

ARAME 

TERS 


















2 

BMIR 

3.15 

15.31 

41.36 

339 

.942 

-0.5 

+18.0 

3.44 

10.59 

6.9 

307 

19 

.76 

1109 

a 

2.25 

1099 

4 

3.38 

356 

16 

.43 

3 

16.3) 

9,45 

14.53 

28. 3D 

132 

.970 

+ 1.4 

-6.6 

5,08 

11.54 

11.6 

165 

19 

.51 

803 

8 

1.57 

825 

4 

2.38 

161 

16 

.21 

6 

LZC 

4.00 

16.67 

41.42 

328 

.954 

+0.3 

+ 168 

4.29 

11.57 

7.4 

254 

19 

.83 

1045 

a 

2.30 

1044 

4 

3.24 

199 

16 

,25_ 

7 

(6.0) 

12.00 

12.59 

28.21 

157 

.958 

0.0 

-7.8 

4,50 

10.64 

8.9 

159 

19 

.45 

814 

8 

1.51 

856 

4 

2 43 

274 

16 

.34 

10 

ETLF 

0.125 

15.42 

34.26 

211 

.959 

+ 1.1 

+6.3 

s.ia 

11.68 

9.7 

225 

19 

.62 

957 

8 

1.89 ' 

960 

4 

—265.— 

420 . , 

— L6-. 

. -JLU. 
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Table 13. Sensitivity Analysis of Little Chestuee (Part 3 of 4) 
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Table 13. Sensitivity Analysis of Little Chestuee (Part 4 of 4) 
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Table 14. Summary of Sensitivity Analysis for Little Chestuee 
(Part 3 of 3) 
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3. 6 CORRELATION (TASK 6) 


The data base for the Correlation task consists of (1) the ten sets 
of ’'best" model parameters from the Calibration task and (2) ten sets of 
Watershed Parameters (i. e. , observable or inferable characteristics) 
from the Physiographic Data Base Task, as indicated in Figure 28. The 
objective of the task was to derive a statistically based set of formulas 
for quantifying those parameters which are normally estimated by the 
calibration process previously described. 

A linear correlation coefficient was calculated for each of 88 
pairs of parameters. In each pair, one parameter was chosen from 
the set of calibrated model parameters and the other from the set of 
"observable” parameters. These correlation coefficients, for a linear 
relationship (y = a + bx), were found to be generally rather low, the 
exceptions (up to 0.8) being associated with parameters whose values 
lie in narrow ranges. The results of this step provided a basis for 
selection of combinations of observable parameters for the multi- 
variate linear regression analysis . It was necessary to accept linear 
correlation coefficients as low as 0.39 in order to generate a useable 
expression for each calibrated parameter. 

1 6 

The result of the multivariate regression analysis 1 is the set 
of equations listed in Table 15. They were used to estimate model 
parameters for the Validation task. The parameters FIMP, OFSS, 
OFMN, OFSL, and OFMNIS can be determined from remotely sensed 
image data, provided resolution (not yet determined precisely) is ade- 
quate and relative elevations can be discerned. The soil parameters 
AWC (Available Water Capacity) and PERM-A (Permeability of the A 
Horizon) were obtained from County Soil Surveys of the U.S. Soil Con- 
servation Service; they are not remotely sensible at present. In a 
future application in an ungaged area, some ground truth data related 
to soil characteristics will be required unless those characteristics 
can be inferred from surrogate information*”’ . 

The data base sample points are widely scattered in all cases. 
This is evident from the statistical indices calculated and shown in 
Table 16. It would obviously be worthwhile to refine the Correlation 
Analysis. Additional observables should be introduced, carefully chosen 
with regard both to hydrologic phenomena and remote sensing system 
capabilities. Additionally, interaction between independent variables 
should be resolved through translation into an intermediate set of com- 
posite variables. Nevertheless, this simplified analysis produced re- 
sults which led to acceptable validation runs and a strong indication of 
the feasibility of the technique. 
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OBSERVABLE AND 
INFERABLE WATERSHED 
CHARACTERISTICS 
FROM TASK 1 


MODEL 
PARAMETERS 
FROM TASK 4 



Figure 28. Task 6: Correlation 
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Table 15. Correlation Results (Multiple Regression Analysis) 


BMIR 

= 

5.27 

4* 

0.75 PERM A 





LZC 

= 

1.80 

+ 

4.73 OFSS 

+ 

2.63 AWC 

+ 

0.26 PERM A 

ETLF 

= 

0.05 

+ 

1.49 OFMN 





SUZC 

= 

0.46 

+ 

2.51 FIMP 

— 

0.04 AWC 



SI AC 

= 

0.80 

+ 

3.84 FIMP 

— 

0.02 AWC 

— 

0.37 OFSS 

BUZC 

= 

0.18 

+ 

9.69 OFMN 

— 

0.0003 OFSL 



BIVF 

= 

-0.02 

— 

9.25 FIMP 

+ 

0.0004 OFSL 



CSRX 

= 

0.95 

+ 

1.92 OFMNIS 

— 

0.003 AWC 



FSRX 

= 

1.01 

- 

0.006 PERM A 

— 

0.00003 OFSL 



BFRC 

= 

0.97 

— 

3.29 FIMP 





BFNLR 

= 

0.86 



4.38 FIMP 






0.0003 OFSL 
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Table 16. Statistical Characteristics, Regression Equations 


PARAMETER 

S.M. 

S.E.E. 

S.M.C.C. 

BMIR 

6.94 

2.89 

0.15 

LZC 

4.96 

1.26 

0.58 

ETLF 

0.135 

0.060 

0.17 

SUZC 

0.290 

0.153 

0.42 

SIAC 

0.300 

0.227 

0.47 

BUZC 

0.349 

0.109 

0.82 

BIVF 

0.296 

0.332 

0.36 

CSRX 

0.957 

0.016 

0.50 

FSRX 

0.959 

0.022 

0.43 

BFRC 

0.935 

0.088 

0.24 

BFNLR 

0.815 

0.097 

0.31 


S. M. = Sample Mean 

S.E.E. = Standard Error of Estimate 

S.M.C.C. = Square of Multiple Correlation Coefficient 



3.7 VALIDATION (TASK 7) 


Five basins of the Tennessee Valley were identified for the 
Validation task. They were exempt from calibration and treated as 
"ungaged" watersheds. The watershed description and historical 
data from the integrated data base, and the model parameters esti- 
mated by the linear regression equations generated in the Validation 
task, are inputs, as indicated in Figure 29. For each basin, a multiyear 
simulation was run, the synthesized streamflow results compared with 
observed streamflow, and analyses of the results printed out in tabular 
and plot formats . 

Results of the validation runs are summarized in Table 17. Mean 
monthly streamflow correlation coefficients indicate generally good agree- 
ment between simulated and observed values. For each basin, the largest 
winter storm event was chosen for comparison of simulated streamflow 
peaks with observed peaks. 

Examples of daily and hourly plots are shown in Figures 30 and 
31 respectively. In the latter, the observed storm event is represented 
only by the flood peak. Peak magnitude and timing are the basis for 
judging the effectiveness of the model with respect to storm events. 
Plotting an observed runoff hydrograph at one-hour intervals would 
require a prohibitively expensive manual translation of strip- chart 
records into hourly runoff for every selected storm event. 

The result of the Validation task is a strong indication of the 
feasibility of making streamflow forecasts for an ungaged basin, at 
least with respect to mean daily streamflow. With respect to storm 
runoff peaks, the result is inconclusive. Further refinement of the 
technique will lead to the ultimate method depicted in Figure 32. 

3.8 RESOLUTION/ ACCURACY ANALYSIS (TASK 8) 

This task was not performed in the study because of budgetary 
limitations. It is being addressed in a separate NASA study, the results 
of which will be available in the second quarter of 1974. 
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Figure 30. Example of Mean Daily Streamflow Plots 
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Figure 31. Example of Hourly Streamflow Plots 
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Figure 32. Ultimate Method 





SECTION 4 


RECOMMENDED ADDITIONAL RESEARCH 


As in all such research, this feasibility study revealed several 
topics deserving further study. Also, the investigators believe that 
some tasks should be done more thoroughly and/ or more efficiently 
with the benefit of experience gained in subsequent tasks. The follow- 
ing paragraphs summarize topics the authors recommend for additional 
research. 

4. 1 EXPANDED CORRELATION DATA BASE 


The five validation runs used data derived from the results of a 
correlation analysis which was in turn based on calibration of only ten 
basins. The size of the sample population limits confidence in the re- 
sults; the good correlations scored in the validation rims may well have 
been n good luck". It is suggested that the following tasks be performed. 

• Collect historical and physiographic data for twenty additional 
watersheds of the Tennessee Valley, including five from the 
Mississippi Alluvial Plain Region. 

• Preprocess and format the new data and integrate it into 
the Master Data Bank. 

• Rework the Correlation Analysis using the expanded data 
base (15 basins). 

• Quantify parameters for ten basins, the five used for vali- 
dation in the basic study and five reserved for the purpose 
from the new set. 

• Perform ten validation runs and analyze the results. 

• Revise this report or issue a supplementary report. 

4.2 RE FINED CORRE LATION STUDY 


Even with the limited population sample (ten sets of parameters) 
used in the basic feasibility study, the correlation analysis should be 
redone with a more rigorous methodology to achieve more refined re- 
sults. It is desirable, for instance, to remove the sub- surface aspects, 
such as permeability and water holding capacity, because they are not 
directly measurable or inferrable from remote sensing. Additional 
physical characteristics which lend themselves to determination by 
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remote sensing, even from satellite altitudes, should be introduced 
into the study.. Hydrologic expertise should be brought to bear in the 
selection of combinations of observable characteristics . Intermediate 
or composite variables should also be derived in order to eliminate the 
effects of interaction among the selected independent parameters. This 
additional research task should produce a report which is supplementary 
to this one . 

4 . 3 VALIDATE TE CHNIQUE IN OTHER REGIONS 


A similar study could be performed in another region which is 
physiographically and climatically very different from the Tennessee 
Valley, as a further validation of the feasibility of the technique. This 
additional work should concentrate in an area in which snowfall is a 
significant form of precipitation. 

4.4 DOCUMENT THE SYSTEM FOR SIMULATION AND ANALYSIS 
OF WATERSHEDS 


The system for simulation and analysis of watersheds, developed 
by IBM around the KWM and OPSET programs, is an effective tool for 
research and hydrologic modeling, river forecasting, and water re- 
sources management. However, it is presently useful only to those 
who are involved in the development of the system, who are familiar 
with its operation. Inadequate documentation exists to transfer the 
tool to another investigating team without an extensive period of learn- 
ing and familiarization. It is suggested that the system for simulation 
and analysis of watersheds be documented by the preparation of operating 
instructions and users' manuals containing flowcharts, program descrip- 
tion, and step by step instructions of how another user could take advantage 
of this system of programs and methods. 

4. 5 PILOT DEMONSTRATION, OPERATIONAL RIVER FORECAST 
APPLICATION 

Accepting the results of the study as a strong indication of the 
feasibility of the technique the study originally set out to investigate, 
it appears appropriate for NASA to leave further refinement of the 
technique to potential users and other researchers. A more appro- 
priate step for NASA to take may well be pursuit of a study jointly 
with a potential user agency, aimed at implementation of a pilot 
demonstration of operational river forecasting utilizing the models, 
methods, and other results of the feasibility study. The first product 
of such an undertaking would be a specification for the demonstration 
system, to be followed by implementation of the pilot system. For a 
period of time the pilot system could be operated in parallel with the 
user agency's normal river forecast methods. Then the two systems 
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could be compared with respect to accuracy and efficiency. Should the 
demonstration prove successful, complete documentation of the system, 
as a NASA product available to potential user agencies, would be the 
final task. 
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APPENDIX A 


SUMMARY OF STUDY HISTORY 


The feasibility study reported herein began in mid- 1971. At the 
beginning of 1972, it was reduced to a low-level sustaining effort and 
then resumed at an appropriate manning level in July 1972. It is 
appropriate to summarize those three segments of the study in this 
Appendix. 

A. 1 LAST HALF, CALENDAR YEAR 1971 

During the months July through December 1971, with the support of 
the Research Institute of the University of Alabama in Huntsville, IBM 

• Selected a family of 55 Tennessee Valley Watersheds, with 
the advice and assistance of the Tennessee Valley Authority 
(TVA) 

• Collected 30 percent of historical data needed for the study 

• Acquired 80 percent of the applicable topographical maps and 
photographic coverage 

• Programmed and operated the statistical rainfall- runoff 
correlation technique sometimes referred to as the "API 
Method" 

• Programmed and operated three basic synthetic hydrograph 
regeneration methods, those of Clark, Nash, and Decoursey, 
with two additional variations on the Nash method 

• After a review of available simulation models, selected the 
Kentucky Watershed Model and its companion calibration 
program for continuation of the study. 

• Established a mutually beneficial technical liaison with the 
Tennessee Valley Authority and orally reported the study 
concept and progress to other organizations: World Bank; 

U. S. Army Corps of Engineers; U. S. Departments of 
Interior (Geological Survey), Agriculture (Agricultural 
Research Service), and Commerce (National Oceanic and 
Atmospheric Administration); and National Academy of 
Sciences. 
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A. 2 FIRST HALF, CALENDAR YEAR 1972 


The initial funding authorization expired at the end of 1971, and the 
study was reduced to a sustaining effort funded at a one-man level from 
January through June 1972, Highlights of this interim effort were as 
follows: 

• The Kentucky Watershed Model and its optimization program 
were acquired and operated in the IBM-Huntsville facility. 
Additionally, the Tennessee Valley Authority Daily Flow 
Model, an experimental simulation program, was acquired, 
set up, and run. 

• Modifications and improvements in the KWM and OPSET 
were begun during this interim period in order to adapt it 
to the purposes of the feasibility study. These additions 
and improvements included preprocessing routines for 
input data and implementation of several plot and tabular 
output options. 

• The historical data base was expanded by the acquisition 
of digital daily streamflow tapes from USGS and digital 
precipitation data tapes from the National Weather Service. 

In order to provide hourly precipitation data for the master 
data bank, a computer program to merge and synchronize 
rainfall data from several hourly and daily precipitation 
gages was completed and made operable. 

A. 3 SECOND HALF, CALENDAR YEAR 1972 

In July 1972 the study was resumed at a five- man funding level, 
which was reduced to three in mid-September. Highlights of the last 
six months of 1972 are as follows: 

• The design of an integrated data bank consisting of historical 
data and model parameters derived from the physiographic 
data base was completed. A master data tape was established 
and data pertaining to three TVA watersheds were entered on 
the tape. The master data tape has the capacity for storage 

of up to 50 watersheds, with data related to 10 years of historical 
experience for each one. 

• The calibration program and simulation model were tested 
using actual data from TVA watersheds. After calibration 
the simulation model is capable of synthesizing streamflows 
which are accurate within five percent of observed flows. 
Accurate simulation of flood peaks with respect to magnitude 
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and time of peak require additional manual adjustments after 
model parameters have been estimated by OPSET. 

• A sensitivity analysis task was introduced into the study in 
order to determine what effects on accuracy of simulation 
variations in the different model parameters have. This 
analysis was performed on two small watersheds of the 
test area. It helped expedite calibration of the remaining 
watersheds. 

• The data base and models were implemented in the IBM facility 
with a Terminal Assistance Package which allows an operator 
to set up programs, modify the data base, and initiate program 
operations from a terminal without having to handle several 
thousands of data cards. 

A. 4 CALENDAR YEAR 1973 

The study was continued at a three-man funding level from January 1 
through September 21, 1973, and was completed, with the results as 
summarized in Section 1. 5 and detailed in Section 3. Most of the effort 
during this time was spent on formatting and integrating the data base for 
15 watersheds. Calibration of ten watersheds required the next largest 
proportion of time, after which the correlation analysis and validation runs 
were quickly completed. 
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